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We measure the sphaleron rate (topological susceptibility) of hot SU(2) gauge theory, using a lattice imple- 
mentation of the hard thermal loop (HTL) effective action. The HTL degrees of freedom are implemented by an 
expansion in spherical harmonics and truncation. Our results for the sphaleron rate agree with the parametric 
prediction of Arnold, Son and Yaffe: F oc a^T^. 



1. MOTIVATION 

Baryon number is not a conserved quantity in 
the Standard Model: due to the anomaly, the vi- 
olation is related to the (Minkowski time) topo- 
logical susceptibility of the SU(2) weak group. 
While at low temperatures the violation is to- 
tally negligible |l[], at temperatures above the 
electroweak symmetry restoration temperature 
(~ 100 GeV) the rate of the baryon number vi- 
olation {sphaleron rate) T is large. This can have 
significant repercussions for baryon number gen- 
eration in the early Universe, and it opens the 
avenue for purely electroweak baryogenesis. 

Even though the weak coupling constant is 
small, at high temperatures the sphaleron pro- 
cesses are dominated by IR momenta k ~ g^T 
and are thus inherently non-perturbative. More- 
over, the IR modes behave essentially classically, 
which is signalled, for example, by the large occu- 
pation numbers of the Bose fields: n{k g^T) — 
[e^l^ - l)-i w T/fc - l/g2 > 1. This has moti- 
vated the much utilized method of using the clas- 
sical equations of motion to calculate P in hot 
SU(2) theories Q (the Higgs and fermionic de- 
grees of freedom efi^ectively decouple in the hot 
EW phase). For recent reviews, see 0, [Q. 

The success of the classical method hinges on 
the efficient decoupling of the almost-classical IR 
modes relevant for the sphaleron processes and 



the strongly non-classical UV modes. However, as 
argued by Arnold, Son and Yaffe this decou- 
pling is not complete. A step beyond the classical 
approximation is the hard thermal loop (HTL) ef- 
fective theory which incorporates the leading 
order effects of the UV modes. The HTL the- 
ory can be cast in various forms; most practical 
for lattice computations is the one where the the 
hard modes are described by including a large 
number of classical massless particles with adjoint 
charge moving on the background of IR fields. 
This field -I- particles system can be put on a lat- 
tice as such, and it has been succesfully used in 
simulations (t). In this work we use an alterna- 
tive Boltzmann- Vlasov approach, where the par- 
ticles are described with local density functions 
n{t,x, k). For full description, see 

2. HTL THEORY ON THE LATTICE 

Let us consider a system consisting of the HTL 
particles moving on the background of IR gauge 
fields. The particle density functions n(t, x, k) 
obey the Vlasov equation 



dt 







daSn + V ■ DSn + dok 



dn 
dk 



DSn 



(1) 



* Presented by K. Rummukainen at the conference LAT- 
TICE '99, Pisa, Italy, July 1999. 



where no — (e*^/-^ — 1) ^, n — no + Sn"", and the 
Lorentz-force (x v x B has been neglected. The 
IR gauge fields evolve according to the Yang-Mills 
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equations: 



Ko = 8.2 gr |3l = 8.7 V=24' 



(2) 



where the 4-velocity v ~ {l,k/k). These equa- 
tions can be further simphfied by factorizing 
(5n° — —gW'^{x,v){dno/dk) and integrating over 
the amphtude |fc| 

v^'D^Wix,v) = v,F°' (3) 

Here d^l integration is over the directions of the 
4-velocity v. The field W'^{x,v) is proportional 
to the flux of the particles at point x to direction 

V. 

In order to perform lattice simulations the field 
W has to be regularized in space (standard lat- 
tice) and on the w-sphere. We do this by ex- 
panding W in spherical harmonics: W°'{x,v) — 
Wi'^{x)Yim{v), and truncating the expansion to 
I ^ ^max- In terms of Wim, the equations (^) 
finally become S 
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Here the coefficients Cl'^l = J driy^^i^'F//™/ and 
~ J dQYimv\ Eq. is the Gauss law, and, 
as long as it is satisfied by the initial configura- 
tion, it is preserved by Eqs. (H) and (^). 

With a finite ^max, these equations can be read- 
ily discretized: SU(2) gauge field is defined on the 
links of the lattice, and the (^max + 1)^ adjoint 
fields are on lattice sites. The discretization 
and the properties of the theory on the lattice are 
discussed in detail in ||]. 

3. THE SPHALERON RATE 

The measurement of the sphaleron rate pro- 
ceeds along similar lines to the purely classi- 
cal theory: first, we generate an ensemble of 
initial thermalized configurations (which satisfy 
the Gauss law), and then evolve these with 
Eqs. (H),(|l)- We then obtain F by measuring the 
rate of the Chern-Simons number diffusion |^,||j3) . 
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Figure 1. The sphaleron rate dependence on the 

^max CUtofi'. 



We have to check how F depends on (a) Zmax, 
(b) lattice spacing and (c) the Debye mass to_d. 
Only the last parameter is physical (it depends 
on the particle content of the theory). 

Let us first consider the ^^ax dependence. In 
Fig. we show F measured from a set of lattices 
with Zmax ^ 10. We note that when Zmax is even, 
the rate remains remarkably constant (much bet- 
ter than indicated by naive arguments Q). The 
behaviour at odd Zmax can be understood by con- 
sidering the properties of the gauge field propa- 
gator 1^. Thus, we conclude that modest values 
of ^max ~ 4-6 are sufficient in order to obtain the 
^max ^ oo behaviour within reasonable statistical 
errors. 

Dimensionally, one would expect that F oc 
a^T^, the non-perturbative scale to the fourth 
power. However, as argued by Arnold, Son and 
Yaffe 1^, the evolution of the IR fields is Landau 
damped by the UV modes, and the rate is slower 
by one further factor of a, parametrically 

T = K'^ct'T^, (7) 



where k' is a constant to be determined by lattice 
measurements. In Fig. ^ we show the behaviour 
of F against g^T'^/m]^, measured using various 
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Figure 2. The sphaleron rate F in HTL theory. 
In leading order /3o = where a is the lattice 

spacing. 

lattice spacings a oc I/Pq. The rate is clearly 
not constant when the Debye mass is varied, and 
it goes to zero when mj~, — > cxi, as predicted by 
Eq. (^). We have not observed any significant 
dependence on the lattice spacing (provided that 
it is small enough). It should be noted that the 
physical mjj is not only the 'bare' mj^ which ap- 
pears in Eq. (|^), but it is a sum of the bare mj^ 
and a contribution due to the UV lattice gauge 
field modes (xl/a |p^ . 

In Fig. H we compare the coefficient k' of the 
scaling law (0) as measured in this work, with the 
particles method 0, and with only the classical 
SU(2) gauge theory evolution without any added 
HTL degrees of freedom |0 . In the last case the 
physical mj^ arises solely through the lattice UV 
modes; here the lattice spacing is up to a factor 
of ~ 4 smaller than in the two HTL approaches. 
The consistency of the results is remarkable, con- 
sidering the very different treatments used. 

To conclude with, the sphaleron rate in hot 
SU(2) gauge theory is now settled. Inserting the 



Standard Model value of 



1 1/6. 9^T^ 



we 



obtain for the rate a value F = (25 ± 2)a^T'^. 
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